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ABSTRACT 



We report the discovery of a new detached eclipsing binary known as 
NSVS01031772 with component masses M x = 0.5428 ± O.OO27M , M 2 = 0.4982 
± 0.0025M o , and radii R x = 0.5260 ± O.OO28i? , R 2 = 0.5088 ± O.OO3Oi? . The 
system has an orbital period of 0.3681414(3) days and an apparent magnitude V 
~ 12.6. The estimated effective temperatures of the stars are T e ffi — 3615 ± 72A" 
and T e ff 2 = 3513 ± 31i^. The space velocities of the system suggest that it is a 
main sequence binary and with a metallicity that is approximately solar. The 
two stars in this binary are located in a region of the Mass-Radius relation where 
no accurate observational data was previously available. Similarly to other low- 
mass binaries recently studied, the radius of each star in NSVS0 1031 772 exceeds 
the best evolutionary model predictions by about 8.5% on average. 

Subject headings: binaries: eclipsing — binaries: spectroscopic — stars: funda- 
mental parameters — stars: late-type — stars: individual (NSVS01031772) 
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1. Introduction 

Low-mass stellar structure models have improved fast in the past few years as reviewed 
by Chabrier & Baraffe (2000) and Chabrier et al. (2005). However, that improvement has 
not yet been matched by observations. Accurate parameters of low-mass stars are difficult to 
obtain, with the best source of precise data being double-lined eclipsing binaries (DDEBs), 
but those systems are scarce. Only two low-mass DDEBs had been found until 1998, CM Dra 
(Lacy 1977; Metcalfe et al. 1996) and YY Gem (Leung & Schneider 1978; Torres & Ribas 

2002) . A few more have been discovered since then, CU Cnc (Delfosse et al. 1999; Ribas 

2003) , BW3 V38 (Maceroni & Montalban 2004), TrES-HerO-07621 (Creevey et al. 2005), GU 
Boo (Lopez-Morales & Ribas 2005), 2MASS J05162281+2607387 (Bayless & Orosz 2006), 
UNSW-TR-2 (Young et al. 2006), and a new binary found in the open cluster NGC1647 by 
Hebb et al. (2006). However, the number of low- mass DDEBs is still small, and the error 
bars on the parameters of some of those systems are too large to place rigorous constraints 
on the evolutionary models. There are also portions of the Mass-Radius relation where no 
good observations exist (e.g. 0.25-0.35 and O.5-O.6M ). 

Clearly, the sample size of well-studied low-mass binaries needs to be increased. Over 
the past two years we have searched for candidate low-mass binaries in the Northern Sky 
Variability Survey (NSVS) database (Wozniak et al. 2004). Our search algorithm uses two 
period-searching methods: the "string/rope" method based on the Lafler-Kinman statistic 
(Clarke 2002), and the analysis of variance method (Schwarzenberg-Czerny 1989). The first 
candidate system, known as NSVS01031772 (hereafter NSVS0103), with J2000 coordinates 
of a — 13:45:35, 5 = +79:23:48 (Figure 1), was identified as a possible low-mass binary 
based on its short orbital period (P m 0.368 days), the duration of its eclipses (~ 16% of 
phase), and its red near-infrared colors, obtained from the 2 Micron All Sky Survey (2MASS) 
database (Skrutskie et al. 2006). In this work we present follow-up spectroscopic and pho- 
tometric observations which confirm the low-mass nature of NSVS0103. We derive accurate 
fundamental parameters for the component stars and compare our results to evolutionary 
models. 



2. Radial Velocity Curves 

We collected a total of 108 spectra in two nights during May 2005 with the echelle 
spectrograph at the 4-m Mayall telescope at the Kitt Peak National Observatory (KPNO). 
The wavelength coverage of each spectrum is 5700-8160 A, with a resolving power of 18,750 
at 6000 A and an average signal-to- noise ratio (S/N) of 7-10 per pixel. We also obtained a 
high S/N spectrum of the M dwarf GJ 740 (M1.5V) for use as a template in the derivation 
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of the radial velocities. By using a real star as template we avoid the problems that the 
low-mass stellar atmosphere models have reproducing some spectral features of the stars. 

The radial velocities were extracted using the implementation of the cross-correlation 
algorithm TODCOR (Zucker & Mazeh 1994) kindly supplied by Guillermo Torres. The 
analysis of NSVS01031 with TODCOR was analogous to those described by Torres & Ribas 
(2002) and Lopez-Morales & Ribas (2005) in their analyses of YY Gem and GU Boo. The 
template was rotationally broadened by 70 km s _1 to optimize the results of TODCOR. 
The standard errors of the resulting radial velocities are 10 km s _1 or less. We tested for 
systematic effects in the derived velocities and found no obvious ones. The radial velocities 
for both stars phased on the eclipse ephemeris derived below are shown in Figure 2. The 
data on that figure are available online. A sample of the contents of the online table is shown 
in Table 1. 



3. Light Curves 

We obtained complete V, R, and /-band light curves of NSVS0103 over 12 nights be- 
tween March and May 2005, using the Apogee U55 512x1152 CCD on the Southeastern 
Association for Research in Astronomy (SARA) 0.9-m telescope at KPNO. The data were 
reduced and analyzed using standard aperture photometry packages in IRAF 1 , with no dif- 
ferential extinction effects taken into account given the relative small separation between 
the target and the comparison and check stars in the field (see Figure 1). The final light 
curves, which contain 649 points in V, 843 points in R, and 1210 points in /, are illustrated 
in Figure 3. The average photometric precision per data point is 0.017 mags in V, 0.013 
mags in R, and 0.007 mags in /. 

NSVS0103 was also observed on the nights of June 29 and 30, 2006 using the 0.6m 
telescope at the Mount Laguna Observatory (MLO), equipped with an SBIG CCD and R 
and / filters. IRAF was again used to apply the flat-fielding and dark current corrections, 
and to derive differential light curves. The MLO light curves, which contain 366 points in 
R and 364 points in /, are shown in Figure 4. The average photometric precision per data 
point is in this case 0.009 mags in R and 0.0011 mags in /. All the data shown in figures 3 
and 4 are available online. Table 2 shows a sample of the contents of the online table. 

Times of minimum light derived from our photometry and from the NSVS database are 



1 IRAF is distributed by the NO AO, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the NSF. 
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given in Table 3. 

From those times of minima we derive the following ephemeris equation 

T MinI = HJD 2, 453, 456.6796(2) + 0.3681414(3)^ (1) 

With the exception of one point, the typical (O—C) residuals are less than one minute. There 
are no obvious trends. Using the orbital period above, we find an average phase difference 
between primary and secondary minima of A0 = 0.4989±0.0027, which is consistent with a 
circular orbit. 



4. Analysis 

Each set of light curves was simultaneously modeled with the velocity curves using 
the ELC code (Orosz & Hauschildt 2000) and its various optimizers using updated model 
atmospheres for low- mass stars and brown dwarfs (Hauschildt, priv. comm.). The model has 
two main types of free parameters: those related to the geometry of the binary such as the 
stellar masses, radii, and separation, and those related to the radiative properties of the stars 
such as their effective temperatures, gravity and limb darkening, and locations of spots (if 
any). For the geometrical parameters, we assumed a fixed orbital period of 0.3681414 days, 
zero eccentricity, and synchronous rotation for both stars. The free geometrical parameters 
were the inclination i, the mass and radius of the primary, Mi and Ri, the ratio of the 
radii R2/R1, and the A-velocity of the primary K\. For a given inclination i and orbital 
period P, the orbital separation a and the mass ratio Q = M 2 /M 1 can be found specifying 
Mi and K x . The value of a sets the absolute scale of the binary, while Q gives a unique 
Roche geometry. We find that Mi and Ki are more efficient to optimize the solutions for 
well-detached binaries than a and Q, since Ki can be inferred from the radial velocity curve 
and Mi, roughly, from the spectral type. Once the scale and the mass ratio are known, 
specifying Ri and the ratio of the radii R2/R1 gives the specific values of the equipotentials. 

For the radiative properties of the stars, we eliminate the need for a parameterized limb 
darkening law by using tabulated model atmosphere with specific intensities over a wide 
range of temperatures and gravities. The gravity darkening exponents were set according 
to the temperatures of the stars following Claret (2001). The stars in NSVS0103 are nearly 
spherical so there is essentially no gravity darkening. We used "simple reflection" (see Wilson 
1990). The free radiative-properties parameters were the temperature of the primary 7\, the 
temperature ratio T 2 /Ti, and the parameters to describe two spots on the primary (this is the 
model that provided the best results in both the SARA and MLO datasets, see below). The 
spots in ELC are specified by a temperature factor, the longitude and latitude of the spot 
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center, and the angular radius of the spot, similarly to the Wilson-Devinney code (1971). 
Finally, ELC has a phase-shift parameter to account for small errors in the ephemeris. In 
total, we have 16 free parameters. 

We began by modelling the radial velocity curves simultaneously with the SARA light 
curves since these light curves have better sampling and more filters than the MLO light 
curves. It became immediately clear that the model radial velocity curves did not match the 
observed radial velocities near the conjunction phases (e.g. during the primary and secondary 
eclipses). During partial eclipse, the spectral line profiles may no longer symmetric, which 
may lead to a difference between the measured radial velocity and the actual radial velocity 
of the star. ELC computes the change in radial velocity (generally known as the Rossiter 
effect) using the technique outlined in Wilson & Sofia (1976), which accounts for the shift 
in the "center of light". Since this simple technique may not fully mimic the way TODCOR 
measures velocities during partial eclipse, we did not include in the fit the primary and 
secondary radial velocities between phases 0.425-0.575 and 0.925-1.075. In all, a total of 38 
points were excluded from each radial velocity curve, leaving 71 points per curve. We note 
that since the mass constraints come mainly from observations near the quadrature phases, 
excluding the points near the conjunction phases will have very little if any effect on our 
results. The points not included in the fit are represented as open symbols in figure 2. 

The SARA light curves are not symmetric about phase 0.5. Usually asymmetries in the 
light curves are attributed to the presence of spots on one or both stars. We tried models 
with a single spot on the primary and models with a single spot on the secondary, but could 
not find an acceptable solution. We then tried models with two spots on the primary, models 
with two spots on the secondary, and models with a spot on each star. The best solution 
was found with a model with two bright spots on the primary. That model is illustrated in 
the top diagram of Figure 5. We run the ELC genetic optimizer for 700 iterations, using the 
final solution to scale the error bars so that \ 2 — N — 1 for each data set. The resulting 
mean errors per point were 0.017 mag in V, 0.013 mag in R, and 0.006 mag in /. After 
scaling the error bars, we ran the genetic code and also a simple grid search optimizer to 
establish the la errors on the fitted and derived parameters using brute-force (see Orosz et 
al. 2002). The results are summarized in the second column of Table 4. Figures 2 and 3 
show the best model radial velocity and light curves plotted with the observations. 

The MLO light curves were modeled in the same way as the SARA light curves described 
above. The results are summarized in the third column of Table 4. Figure 4 shows the best 
model fits for those observations. The two bottom diagrams in figure 5 show the best spot 
configuration in this case. The light curves from MLO turned out to be not as useful as the 
SARA light curves for deriving precise parameters of the system. The statistical errors on 
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the fitted and derived parameters from the MLO data are typically a factor of three larger 
than the statistical errors from the SARA data. The SARA light curves have a total of 2702 
points in three filters obtained with a 0.9m telescope, whereas the MLO light curves have a 
total of 730 points in two filters obtained with a 0.6m telescope, so it is not surprising that 
much more precise results were obtained from the SARA data. 

5. Parameters of NSVS01031772 

We use the ELC solutions discussed in §4 to derive the physical parameters of NSVS0103. 
The masses of the stars derived from the orbital solution of the SARA data are Mi— 0.5416 
± O.OO68M and M 2 = 0.4988 ± O.OO48M . For the radii we obtain R l= 0.5273 ± 0.0029i? o 
and i?2= 0.5058 ± O.OO32i? . The values of those parameters from the MLO data are Mi— 
0.543 ± O.OO3M and M 2 = 0.498 ± O.OO3M and R x = 0.510 ± 0.010i? Q and R 2 = 0.538 
± 0.010i? o . Given the discrepancies between the two datasets, mainly in the values of the 
stellar radii, we decided to adopt a weighted average of those parameters as our final values. 
The results are shown in Table 5. 

The projected rotational velocities of the stars, derived from their radii and the orbital 
period of the system (assuming synchronous rotation), are f synci sini = 72.13 ± 0.38 and 
^sync 2 srn ^ — 69.77 ± 0.41 km s -1 , consistent with the results from TODCOR. The values of 
T e ffi and T e ff 2 in Table 5 correspond to the average of the SARA and MLO results. Note 
that the effective temperatures reported in Table 4 correspond to the temperature of the 
photosphere of the stars without spots. We recomputed the values of T e //i to account for 
the presence of spots. The resultant mean effective temperature of the system is therefore 
T e ff = 3564i73.fr. This value agrees within the errors with the mean effective temperature 
of T e ff = 3750 ± 185K derived using the color-T e // relations in Table 2 of Lopez- Morales & 
Ribas (2005). 

The luminosity and absolute magnitude of the stars were computed from their effective 
temperatures, radii, and surface gravities using filter-integrated NextGen model atmo- 
spheres provided by France Allard. A simple computer code was written that uses the 
temperature, gravity, and radius plus the la errors of each component as inputs to compute 
the absolute magnitudes of the binary in the standard UBVRIJHK filters. We find, for 
example, M v = 9.08 ± 0.22, Mj = 5.77 ± 0.09, and M K = 4.89 ± 0.08. The distance can 
be computed by using the observed near infrared magnitudes given in the 2MASS database, 
where NSVS0103 has the designation 2MASS J13453489+7923482. The observed magni- 
tudes are listed as J = 9.692 ± 0.021, H = 9.021 ± 0.018, and K = 8.778 ± 0.016. The 
.fT-band extinction to infinity in this direction is 0.012 mag, according to the extinction 
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maps (see http://irsa.ipac.caltech.edu/applications/DUST/). Assuming the temperatures 
above, and the best-fitting radii and gravities, we find a distance of d = 60.3 ± 1.0 pc using 
the observed K magnitude and extinction and a distance of d — 60.9 ±0.9 using the J-band 
magnitude and extinction. 



5.1. Age, Space Velocities and Stellar Activity Indicators 

NSVS0103 does not seem to be related to any known cluster, stellar association, or 
star formation region. Therefore, we can only evaluate its age from its space velocities 
(U, V, W) 2 . The heliocentric space velocities of NSVS0103 were computed from its position, 
radial velocity (7 = 19.0 ± 1.0 km s _1 ), distance (d = 60.6 ± 0.7 pc), and proper motions; 
the latter retrieved from the USNO-B1.0 catalog (Monet et al. 2003): fi a = 100 ± 4 and 
/i«5 = 62 ± 1 mas yr _1 . The obtained components of the space motion are U = —21.07 ±0.16, 
V = 16.54 ± 0.04 and W = —1.19 ± 0.15 km s^ 1 , which correspond to a total space velocity 
of S = 26.81 ± 0.13 km s _1 . The value of W indicates that NSVS0103 is confined to the 
galactic plane. The location of the binary in the U — V plane does not fall within the young 
disk stars area defined by Eggen (1989), nor the area occupied by the young population 
tracers (Skuljan et al. 1999). Also the space velocities of NSVS0103 do not match any 
known moving group (Montes et al. 2001). We conclude therefore that NSVS0103 is not a 
young object and has most likely already reached the main sequence. In addition we expect 
the metal abundance of NSVS0103 to be close to solar, since its space motions agree with a 
disk population. 

NSVS0103 appears in the ROSAT All-Sky Bright Source Catalog (Voges et al. 1999) 
as X-ray source 1RXS J134540.6±792332. Using the calibration equation by Schmitt et 
al. (1995) we estimate an X-ray flux of (5.97±1.25) x 10~ 13 ergs cm -2 s -1 , which results in an 
X-ray luminosity of logLx(ergs s _1 ) = 29.39±0.09, which is similar to the X-ray luminosities 
of YY Gem, logL x (ergs s' 1 ) = 29.27 ± 0.02, and GU Boo, logL x (ergs s" 1 ) = 29.3 ± 0.2, 
computed using the same calibration equation. Finally, we observe strong H a emission lines 
in our spectra, with an average equivalent width of 5.4 ±1.1 A over all phases. That H a 
emission level is higher than those of YY Gem (2.0 A), CU Cnc (3.85-4.05 A) and GU Boo 
(1.7 A). 



2 Positivc values of U, V, and W correspond to velocities toward the Galactic center, Galactic rotation 
and North Galactic pole. 
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6. Comparison with Models 

We compare the masses and radii of the stars in NSVS0103 to the predictions by the 
models of Baraffe et al. (1998) (hereafter B98) and Siess et al. (2000) (hereafter S00), which 
are the only models that explicitly attempt to reproduce the properties of main sequence low- 
mass stars. We also include in this comparison the other known binaries compiled in Table 
10 of Lopez- Morales & Ribas (2005), and the new systems found by Hebb et al. (2006), 
Bayless & Orosz (2006) and Young et al. (2006). After exploring a range of metallicities and 
ages, we conclude that the Z=0.02 models provide the closest values to the measured masses 
and radii of NSVS0103. Different isochrones give practically identical results, assuming the 
system has already reached the main sequence 3 , so we chose a 0.35 Gyr isochrone which 
also allows us to compare NSVS0103 to the other binaries. This comparison is illustrated in 
Figure 6. We find that the models consistently predict smaller radii than the ones observed. 
The discrepancies are of the order of 5% for the primary and 11.9% for the secondary for 
B98. In the case of the S00 models, those discrepancies are 15.5% and 22%. This result 
fully agrees with the trend observed in previous studies of YY Gem, CU Cnc, and GU Boo. 
We also show in figure 6 the empirical Mass-Radius relation derived by Bayless & Orosz 
(2006), using the parameters of previously studied low-mass binaries. The parameters of the 
secondary in NSVS0103 fit well in this case, but the radius of the primary is overestimated 
by about 5% 

We compare next the T e gf and My of the stars in NSVS0103 to the two models above. 
The error bars of these parameters are much larger than in M and R, since we have to use 
external calibrations to compute their values. However, they are still useful to check for 
significant discrepancies between models and observations. The temperatures predicted by 
B98 for 0.50 and 0.54 M Q stars are 3657K and 3750K, while the predictions by S00 for those 
same stars are 3810K and 3880K. In both cases the models predict hotter temperatures than 
the ones obtained for NSVS0103. Finally, B98 predict absolute magnitudes of 9.89 and 9.48 
for 0.50 and 0.54 M Q stars, while the values predicted by S00 are 9.97 and 9.56. In both 
cases the absolute magnitudes of NSVS0103 agree with the models, given the uncertainties. 



3 Chabrier & Baraffe (1995) estimate a time of arrival to the ZAMS of about 0.3 Gyr for stars below 0.6 
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7. Summary and Discussion 

We have found a new detached eclipsing binary composed of two M-type stars with 
masses between 0.49 and 0.53 M Q . Both stars fall in a stretch of the mass-radius relation 
where no accurate data were available before. The radii of these stars are the most accurate 
derive to date for low-mass stars, with uncertainties of only 0.55 and 0.63%. The estimated 
temperatures of the stars are T cSl = 3615 ± 72K and T cS , 2 = 3513 ± 31K, after accounting 
for the effect of spots. We also conclude from the space velocities that the binary is most 
likely on the main sequence stage. Therefore the parameters of the stars are suitable to test 
stellar structure models. 

The models that most closely reproduce the radius of the stars in NSVS0103 are B98 
(Z = 0.02). However, they still underestimate the observed values by 8.5% on average; a 
result similar to previous studies of other binaries. This discrepancy between models and 
observations could be attributed to innacuracies in the equation of state. However, the same 
models describe well stars in long period (> 10 days) spectroscopic binaries and single field 
stars (Delfosse et al. 2000; Ribas 2005). One difference between short period binaries and 
the long period binaries and single stars is the higher rotational velocities of the former 
as result of orbital synchronization. The higher rotational speeds enhance the magnetic 
activity, and hence the appearance of large spots, X-ray activity, and emission lines in the 
stars. The observed larger radii and lower temperatures could be a consequence of that 
enhanced activity. Mullan & MacDonald (2001), for example, suggest that the larger radii 
may be caused by the inhibition of convection in stars with strong magnetic fields. The 
stars would appear cooler as a consequence of the larger radii. If this hypothesis is correct, 
it affects not only the components of short period binaries, but also any active low-mass 
star. If the radius of low-mass stars is correlated to their magnetic activity (equivalently to 
their rotational velocity), the current models may only apply to slow, inactive rotators. For 
younger, faster rotating stars, their magnetic activity level may have to be included in the 
models as an additional parameter to compute their fundamental properties. 
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Table 1: Sample of radial velocities table. A full version of this table is available online 



HJD (days) 


phase 


Kx (km/s) 


K\ err (km/s) 


K 2 (km/s) 


K 2 err (km/s) 


2453493.7683 


0.746 


160.51 


9.29 


-134.44 


14.74 


2453493.7710 


0.753 


158.48 


11.32 


-136.25 


9.47 


2453493.7737 


0.760 


164.10 


14.24 


-135.40 


12.21 



Table 2: Sample of SARA and MLO light curve data. A full version of this table is available 
online 



HJD (days) 


phase 


mag 


mag err 


2453456.92852 


0.676 


0.466 


0.018 


2453456.92971 


0.679 


0.458 


0.016 


2453456.93091 


0.683 


0.470 


0.016 


2453456.93212 


0.686 


0.481 


0.016 


2453456.93331 


0.689 


0.464 


0.016 



Table 3: Observed Times of Minimum Light 



Time of minima 


cycle # 


O-C 


(HJD) 


(since To) 


(min) 


2453473.98266 


47.0 


0.596448 


2453478.95221 


60.5 


0.079632 


2453480.97698 


66.0 


0.068544 


2453480.79303 


65.5 


0.242352 


2453490.91543 


93.0 


-1.901088 


2453492.94170 


98.5 


0.247824 


2453493.86213 


101.0 


0.357984 


2453520.73664 


174.0 


0.628416 


2453520.92006 


174.5 


-0.308592 


2453916.67214 


1249.5 


-0.200592 


2453916.85669 


1250.0 


0.489600 


2453917.77661 


1252.5 


-0.134640 
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Table 4: Model Parameters for Light Curve Solution 



Parameter 


Value from SARA 


Value from MLO 


orbital period P (days) 


fixed at 0.368141 


eccentricity e 


fixed at 0.0 


rotational velocities 


fixed at synchronous 


0'ravitv darkening nrimarv 


0.0450244 




0.0450994 


limb darkening coefficients 


from model atmospheres 


reflection 


"simple" 


reflection 


inclination 7 fdet?^ 


85.91 ±0.03 


85.86±0.05 


mass ratio O — M^lM-\ 


0.9217±0.0048 


0.9166 ±0.0050 


radius ratio /?o / /?i 


0.9638 ±0.0066 


1.054±0.051 


temperature ratio T2/T1 


1.01129 ±0.00036 


0.9918±0.0012 


primary temperature T\ (K) 


3505 ±6 


3512±14 


secondary temperature T2 (K) 


3545 ±6 


3482 ±18 


senaration ( 


2.1870 ±0.0089 


2.1910±0.0050 


(km s" 1 ) 


143.85 ±0.37 


143.5±0.6 


A :; (km s" 1 ) 


156.06±0.88 


156.6±0.7 


Omega-potential 1 


5.121±0.015 


5.24±0.11 


Omega-potential 2 


5.053 ±0.022 


4.81 ±0.08 


temperature factor, spot 1 ... 


1.048 ±0.004 


1.19±0.02 


latitude spot 1 (deg) 


26.1±0.3 


104±17 


longitude spot 1 (deg) 


310.9±0.5 


291.5±2.0 


angular radius spot 1 (deg) . . . 


88.4±0.5 


8.5±1.0 


temperature factor, spot 2 ... 


1.196±0.002 


1.075±0.010 


latitude spot 2 (deg) 


42.4±1.0 


46 ±6 


longitude spot 2 (deg) 


90.6 ±0.5 


63.5±2.0 


angular radius spot 2 (deg) . . . 


14.9 ±0.2 


24.5±0.5 



Table 5: Absolute dimensions and main physical parameters of the components of NSVS0103 



Parameter Primary Secondary 

Mass (Mo) 0.5428 ± 0.0027 0.4982 ± 0.0025 

Radius (i? ) 0.5260 ± 0.0028 0.5088 ± 0.0030 

log 5 (cgs) 4.730 ±0.005 4.722 ± 0.006 

v sync sini (km s" 1 ) 72.13 ±0.38 69.77 ±0.41 

T eff (K) 3615 ± 72 3513 ± 31 

L/Lq 0.0426 ± 0.0034 0.0356 ± 0.0013 

M bo i (mag) 8.08 ± 0.25 8.27 ± 0.12 

M v (mag) 9.74 ± 0.22 10.04 ± 0.10 



i : i iJL ' - ! ' i : LMi. i*: r. .-. V " i f I t. t I i ■■ { A, _£j L v ' -^Pt ; 

13 h 47 m 30 s 13 h 46 m 30 8 13 h 45 m 30 8 1 3 h 44 ffl 30 B 13 h « ffl 30 

Right Ascension (J2000) 



Fig. I.— Finding chart for NSVS0103, made from the Digitized Sky Survey. NSVS0103 
marked with the hash marks. 
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Fig. 2. — Radial velocity curves of NSVS0103. The filled circles and squares correspond, 
respectively, to the velocities of the primary and the secondary used in the orbital fit in §4. 
The open symbols show points excluded from that fit. The solid and dotted lines represent 
the orbital solution obtained with ELC. The dashed line shows the velocity of the center of 
mass of the system (~ 19 km/s). These data are available in the electronic edition of this 
journal (Table 1). 
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0.5 1 

Orbital phase 



Fig. 3. — V, R, and I light curves from SARA. The dots represent individual observations in 
each passband. The solid lines show the best fit to the data using ELC. The O-C diagrams 
on the top show the residuals of those fits. These light curves are available in the electronic 
edition of this journal (Table 2). 
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0.5 1 

Orbital phase 

Fig. 4. — R and I light curves from MLO. The dots represent individual observations in each 
passband. The solid lines show the best fit to the data using ELC. The O-C diagrams on the 
top show the residuals of those fits. These light curves are available in the electronic edition 
of this journal (Table 2). 
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Fig. 5. — Spot models that best reproduce the SARA light curves (top diagram) and the 
MLO light curves (two bottom diagrams). In both cases the best solution correspond to two 
bright spots on the primary. The parameters of the spots are listed in Table 4 
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Fig. 6. — Mass-radius relations of stars between 0.1 and 0.8 M & predicted by the models 
of Baraffe et al. (1998) (solid line) and Siess et al.(2000) (dashed line), for an age of 0.35 
Gyrs and Z = 0.02, and the empirical relation of Bayless & Orosz (2006) ((dotted line). 
The open triangles show the location of all the low-mass main sequence binaries previously 
found. The open circles show the two stars in NSVS0103. Each point includes error bars. 



